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We investigate the effect of mass anisotropy on the Wigner crystallization transition in a two- 
dimensional (2D) electron gas. The static and dynamical properties of a 2D Wigner crystal have 
been calculated for arbitrary 2D Bravais lattices in the presence of anisotropic mass, as may be 
obtainable in Si MOSFETs with (110) surface. By studying the stability of all possible lattices, we 
find significant change in the crystal structure and melting density of the electron lattice with the 
lowest ground state energy. 



One of the most unexpected discoveries in the 2D elec- 
tron systems in semiconductor structures is the possibil- 
ity of a metal-insulator transition (MIT) first suggested 
by experiments of Kravchenko et at. in high mobil- 
ity Si metal-oxide-semiconductor field-effect transistors 
(MOSFETs) in zero magnetic field [Q. Later, a num- 
ber of groups have reported possible MITs in 2D elec- 
tron or hole systems in several different semiconductor 
heterostructures, as well as in MOSFETs ||-||. How- 
ever, a genuine MIT in two dimensions is in contradic- 
tion to the scaling theory of localization, which predicts 
that in the absence of electron-electron (or hole-hole) in- 
teractions no true metallic behavior is possible in two 
dimensions with pure potential scattering [[)). In fact, no 
fundamental principle requires that such a scaling argu- 
ment holds in the presence of strong interactions, which 
can be measured by the ratio of the Coulomb energy to 
Fermi energy, expressed by the dimensionless parame- 
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of suggestions that this is the case [ 10 
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cal models not involving MIT have also been proposed 
[]l~3| — . Of all the experimental systems, the critical r s 
varies from 5 in Sio.8sGeo.12 @ up to 35 in GaAs/AlGaAs 
H , which clearly suggests that the Coulomb interactions 
are certainly not negligible, but in fact may be the dom- 
inating scale. 

It is well known that 2D electrons crystallize into a 
triangular lattice (Wigner crystalization) in the low den- 
sity limit where electron-electron interactions dominate. 
In an ideally clean 2D system, the critical r s is predicted 
by Tanatar and Ceperley to be 37 ± 5 from quantum 
Monte Carlo simulations Jl^]. Chui and Tanatar further 
found that the effect of impurities can lower the critical 
r s at which solidification occurs (at least on short length 
scale) to as low as 7.5 fl7|| . This is in good agreement 
with the range of the critical r s observed in the recent 
experiments. 

In this paper, we have carried out a further test on 
the relevance of the Wigner crystallization phenomenon 
to the putative 2D MITs by exploring the effect of 
mass anisotropy on Wigner crystallization, since mass 
anisotropy can be presented in Si- or Ge-based semi- 
conductor devices. In particular, we have studied the 
ground state energy of arbitrary 2D Bravais lattice with 
anisotropic mass. To the best of our knowledge, no result 



on such properties of 2D Wigner crystal with anisotropic 
mass has been published so far. 

The ground state energy of a lattice in the low density 
limit, where exchange process are negligible, can be writ- 
ten as the sum of static Coulomb energy and vibrational 
zero-point energy. With the aid of Ewald's transforma- 
tion [|l9| , we can write the static ground state energy per 
electron of a Bravais lattice as a sum over its lattice sites 
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where <j> n (x) is the Misra function 
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and v is the area of the primitive unit cell of the direct 
lattice. We assume the existence of a neutralizing back- 
ground of positive charges to avoid the divergence in the 
evaluation of the static energy, which is independent of 
the mass anisotropy. In a harmonic approximation, the 
spectrum of lattice vibrations is determined by 
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where u(r°) = r — r° is the deviation from equilibrium 
of the electron whose equilibrium site is r°. The atomic 



force constants <J> a?/ #(r9, r*-) are defined by 
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with v(r) = e 2 /r. The normal mode frequencies w\(q) 
are solutions of the 2D eigenvalue problem: 



M^(q)e A = D(q)e A , 
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where the dynamical matrix D(q) is given by 

D(q)=^$(R )e- 2q - R °. (6) 

R° 

In anisotropic systems, M is a 2 x 2 matrix. The dynam- 
ical ground state energy can be obtained by integrating 
zero-point energies of all modes within the Brillouin zone. 

Bonsall and Maradudin did a comprehensive study on 
the ground state energies of five different crystal struc- 
tures, including triangular lattice and square lattice, for 
isotropic 2D electrons [|l8| . They found that the triangu- 
lar lattice has the lowest energy in low density limit, 
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in units of Ryd = m*e A /h 2 . In Eq. [?], the first term 
comes from the Coulomb interactions between electrons 
sitting on the lattice sites, while the second term is from 
the harmonic oscillations of electrons around their equi- 
librium positions. They also pointed out that the trans- 
verse branch of the dispersion relation is pure imaginary 
for certain directions in the square lattice, implying a 
dynamical instability of this lattice. 
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FIG. 1. 2D Bravais lattices can be mapped to points 
within the shaded area. Dynamically stable lattices are 
mapped to only one tenth of the irreducible area. 

In the presence of the mass anisotropy, the ground 
state may not necessarily be one of the most symmet- 
ric lattices. Therefore, we have developed a systematic 
method to explore the ground state energies of 2D Bra- 
vais lattices with different electron concentrations and 
mass anisotropy. First, we map an arbitrary Bravias lat- 
tice to a point on a 2D plane by the following scheme. 
From an arbitrary lattice site, which is chosen to be ori- 
gin, one can choose an arbitrary pair of two non-collinear 
lattice vectors that span the Bravais lattice. If one ro- 
tates and scales the lattice so that one of the two prim- 
itive vectors lies along x-axis normalized to unit length, 



the end point of the other lattice vector represents the 
2D lattice structure. However, a lattice can be mapped 
to infinite number of points due to arbitrary choice of 
primitive vectors. One can reduce the area of points by 
choosing two shortest vectors from the origin with the 
second shortest one placed along x-axis. Thus, different 
lattices are represented by points in the upper positive 
quadrant as shown in Fig. [I. Using reflection symmetry, 
however, all the points can be confined by the y-axis and 
the two circles shown (the shaded area) in Fig. |l|. The 

three corners (0,1), and (0,0) represent square 

lattice, triangular lattice, and quasi-one-dimensional lat- 
tice, respectively. We can then sample 2D lattices by 
applying a rectangular mesh on the area and calculating 
the ground state energy on each grid point, which can be 
uniquely mapped back to a 2D Bravais lattice. 

Applying Bonsall and Maradudin's calculations on all 
sampled lattices with isotropic mass, we found that most 
lattices have imaginary vibrational modes, therefore are 
unstable. Under our mapping, only lattices around the 
triangular lattice are stable, which occupies roughly 10% 
of the reduced zone, as shown in Fig. [l| Using symmetry, 
one can actually show the triangular lattice is the center 
of all stable lattices. Within the stable area the triangular 
lattice has the lowest energy for r s > 25, which is slightly 
below the Wigner crystallization density r s = 37 ± 5 pre- 
dicted in quantum Monte Carlo calculations |1J. When 
r s = 25, the lattice represented by point A (0.27, 0.95) 
starts to have lower energy than the triangular lattice, 
whose ground state energy is found to be 
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Since point A is adjacent to the unstable area, this im- 
plies that Wigner lattice is no longer the ground state of 
the 2D electron system. 

We have studied different ratios between longitudinal 
mass mi and transverse mass m t . In this paper, we 
present results of mi/mt = 3, which is approximately the 
mass ratio in silicon (110) surface structures We use 
the geometric mean of the two mass components as the 
effective mass in calculating r s . The relative orientation 
of the lattice vectors with respect to the principle axes 
of anisotropic mass has also been taken into considera- 
tion. Figure || compares the dispersion relation of the 
anisotropic triangular lattice with that of the isotropic 
lattice. Mass anisotropy lifts the 2-fold degeneracy at J- 
point thus the longitudinal branch as well. The ground 
state energy for triangular lattice fo mi /m t = 3 therefore 
becomes 
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Compared with Eq. (1), one concludes the mass 
anisotropy leads to an increase in the total lattice vi- 
brational energy. 
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FIG. 2. Dispersion curves for triangular lattice with 
anisotropic mass ratio mi/mt = 3. Frequency is measured 
in units of (ra*e 4 /fc 3 r?^ 2 ). 



density for ideally clean 2D anisotropic electrons, such as 
silicon (110) surface electrons, is much lower than that 
predicted for an isotropic system. 
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FIG. 3. The ground state move away from the trian- 
gular lattice with decreasing r s in anisotropic system with 
mi/mt — 3. 



Since the electrostatic ground state energy is indepen- 
dent of mass anisotropy, one expects that at low enough 
density, the triangular lattice remains as the lowest en- 
ergy configuration. However, just for mi/m t = 3 even at 
very low density, specially r s < 1000, we find that the 
lattice with the lowest ground state energy starts to shift 
along one of the circular curve that confines the reduced 
lattice area towards the quasi-one-dimensional lattice. It 
is understandable that a lattice can stretch along the axis 
of the smaller mass to reduce the longitudinal energy of 
lattice vibration. This trend is shown in the enlarged 
stable area in Fig. ||. 

Fig. |] shows the plots of the ground state energy for 
different lattices represented by the two circular curves 
(A-T-B) at different r s for mi/m t = 3.0. When r s = 100, 
the ground state energy is minimized at a point along 
curve TB. We expect 2D electrons crystallize in the 
lattice structure mapped to this point at such density. 
When r s = 80 the lattice with the lowest energy moves to 
the corner of the stable area, B, which connects with the 
area of unstable lattices. The comparison of the phonon 
dispersion curves of the triangular lattice and the lattice 
mapped to point B suggests that, by crystallizing in less 
symmetric lattice, the electrons gain larger energy in lon- 
gitudinal modes than the loss in transverse modes. When 
r s = 60, the ground state energy decreases linearly when 
approaching point B from the stable area along the cir- 
cular curve TB, which implies that no stable lattice can 
be formed at this density. The comparison of the ground 
state energies suggests that the Wigner crystallization 
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FIG. 4. The ground state energy, in units of Ryd, along 
the circular boundaries (A-T-B) for r s =60, 80, and 100. 

In this calculation we have neglected quantum me- 
chanics of the electrons (i.e. electronic exchange en- 
ergy). While this may have a quantitative effect, it is un- 
likely to completely alter our basic conclusion that r s of 
the Wigner crystallization transition increases with mass 
anisotropy. This is especially true since the r s of interest 
(~ 80) is very large, where exchange energies are small- 
est. 

To summarize, we have found the ground state of 
2D electron system has strong dependence on the mass 
anisotropy. In a pure system, anisotropic electrons re- 
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quire larger spacing to form Wigner lattice so as to sta- 
blize the long wavelength transverse modes. It would be 
interesting to see if this dramatic reduction in Wigner 
crystal density can be observed in extremely clean Si 
(110) MOSFET, in which our calculations estimate the 
critical density, in terms of r s , can be as large as 80. 
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